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Abstract 



A model consisting of a single massless scalar field with a topological cou- 
pling to a pure gauge field is defined and studied. It possesses an SL{2, 
symmetry as a consequence of the gauge invariance. We propose that by 
adding impurities the model can be used to describe transitions between 
Quantum Hall plateaux. This leads to a correlation length exponent of 20/9, 
in excellent agreement with the most recent experimental measurements. 
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I. INTRODUCTION 



Theoretical studies of the Quantum Hall effect can roughly be classified either as studies 
of the many-body physics on a plateau (Laughlin wave function, etc.) or as studies of the 
plateau to plateau transitions. For a review of the latter see [||]. The transitions have been 
viewed as metal-insulator quantum phase transitions and viewed as such are problems in 
Anderson localization |^. The theoretical framework developed to study these problems 
involves the study of electrons in a disordered potential. Averaging over disorder using 
replicas leads to a sigma model in 2 dimensions @]. 

In applying this framework to the Quantum Hall transitions important progress was 
made by Pruisken ||^ 0] who understood that a topological term proportional to the Hall 
conductivity a^y was essential for obtaining the correct features of the phase diagram. For 
the super symmetric version see 0. The topological term modifies significantly the renormal- 
ization group and hence the infra-red properties, in the same way that the 0(3) non-linear 
sigma model with topological coupling at 6' = tt has a non-trivial infra-red fixed point P). 

Despite this progress, the detailed critical properties, such as the exponents, have so far 
remained uncomputable. This is due to the complexity of the sigma models obtained and 
the usual difficulty in finding non-trivial fixed points in the infra-red. It remains unclear 
whether the correct critical exponents are even contained in the sigma-models. 

Progress has been made by inventing a network model which is believed to be in the same 
universality class as the transition |p. Here the transition resembles percolation through 
the impurities. The network model was recently shown ||10| to be equivalent to fermions 
with Dirac-like hamiltonian with various random potentials of the kind studied in [|ll|. For 
these models one is again unable to compute the critical exponents analytically due to the 
complexities introduced by disorder averaging and one has to resort to numerical methods. 

On the other hand, very general and elegant arguments, based largely on gauge invari- 
ance, were given by Laughlin and Halperin to explain the very existence of the plateaux |12 



13| . The argument does not depend on any detailed properties of the disorder; it just has 



to be there to localize some states. Indeed, Halperin argued that from the gauge invariance 
alone, one can infer the existence of both localized and extended states; the argument does 
not involve finding a critical strength of the disorder for a localization/delocalization phase 
transition. What helps the situation is the fact that in the scaling approach to Anderson 
localization, two dimensions is distinguished. In principle, when there is no magnetic field, 
states in two dimensions are always localized, no matter what the strength of the disor- 
der |jl5|. The most extreme conclusion to be reached based on these observations is that 



perhaps even the critical exponents governing the transitions between plateaux may follow 
largely from gauge invariance. It would indeed be very satisfying if both the existence of the 
plateaux and the critical properties of the transitions between plateaux were consequences 
of the same fundamental principle of gauge invariance. In this paper I will construct a model 
based on this notion. In other words, we will simply assume disorder is irrelevant. The issue 
of whether disorder is relevant or irrelevant in the renormalization group will be addressed 
in a separate publication . 

In the approach developed here, in a sense I start from the end (the critical theory) 
rather than from the beginning. Namely, we construct a conformal theory that contains the 
essential features of the extended states. We then add impurities as a perturbation. What is 
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normally considered in the literature is the reverse, i.e. one starts from the localized states 
and looks for a delocalization transition. Critical theories in two dimensions are strongly 



constrained by conformal invariance [|I6|, and the computation of the critical exponents 
should amount to the proper exercise in conformal field theory. This involves identifying the 
appropriate class of conformal field theory. Guided by the importance of gauge invariance, 
the significance of topological terms (the topological term in the sigma-model approach 
should survive in the infra-red since it is proportional to axy), we construct a certain model 
of conformal field theory which possesses the necessary ingredients. The hamiltonian corre- 
sponds to 2-component Dirac fermions. The definition of the model includes a constraint on 
the zero modes that incorporates an essential feature of the 2 + 1 dimensional Quantum Hall 
dynamics. We view this model as describing the essential features of the extended states 
in a pure Quantum Hall system. After bosonizing the fermions and taking the gauge field 
to be a pure, but singular, gauge, we obtain a model with a topological coupling similar to 
the multi-boson models studied originally in the context of string theory [O [ITH]. However 



there are some important differences. Whereas the simplest (smallest) model studied in |T7 
||T8| has Virasoro central charge c = 2 (two scalar fields), the gauge invariance of our model 
effectively reduces the degrees of freedom to a c = 1 theory. This is a new class of conformal 
field theory characterized by an electric/magnetic SL{2,'Z) duality of the kind found by 
Cardy which acts on the modular parameter t = 9/271 + ig'^ where 9 are couplings 



in our model. A number of works have postulated a role for SL{2, Z) in the Quantum Hall 
transitions, based on the phenomenology of transition selection rules [|20| |21|] [^] p4| . 



Though our microscopic model does not have some of the features assumed in these works, 
it may perhaps still be useful towards developing these ideas. Our model is easily studied 
without utilizing the full machinery of conformal field theory; the elementary tools we need 
are bosonization E^ and the properties of exponentials of massless scalar fields, reviewed 



in mn 

Starting from the Kubo formula in 2 + 1 dimensions we derive some formulas for the 
conductivities a^x, o^xy which involve current-current correlation functions in the 2 dimen- 
sional quantum field theory. These formulas are different from the kind of formulas found 
in the literature, the latter being expressed as a double product of retarded and advanced 
1-particle Green functions. We argue that our formulas represent the conformal contribution 
to the conductivity. In our model the conductivities are functions of g, 9, and when axx = 0, 

9 = l/(Jxy 

We then add an impurity potential as a perturbation away from the critical point. If, as is 
normally done, one takes a gaussian random potential, then according to the usual approach 
to Anderson localization, one must average over the disorder and search for a critical point 
corresponding to the critical strength of the disorder for a localization/delocalization phase 
transition. For the reasons given above, we chose the potential so that no disorder averaging 
is needed. In real experiments the critical exponents are measured for a single sample, 
whereas disorder averaging amounts to averaging over different realizations of disorder. If 
we assume that disorder is irrelevant, then universality implies that any impurity potential 
which breaks translational invariance should suffice. Hence, viewing the sample as a very 



^For a more comprehensive review see [27]. 
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large disk, we take the impurity potential to be a circle of impurities somewhere inside the 
disk. This Corbino disk geometry has essentially all of the ingredients used in the gauge 
argument of Halperin [|13|. The virtue of doing this is that the critical properties do not 
depend on the strength of the impurity potential. Furthermore we are already at a fixed 
point before perturbation. Our resulting model can be solved by mapping it onto a boundary 
field theory, with the impurities residing at the boundary. 

In the scaling theory developed in the primary critical exponents are /i, v defined as 



AB oc T^, 



'dp. 



xy 



dB 



oc T' 



and also 



oc |5 - B, 



:i.2) 



Here, pxy is the resistivity, B the magnetic field, Ai? the width of the region between 
plateaux, B^. the critical value at a transition, T the temperature and the correlation 
length. The exponents are related by /i = p/2z/, where p is the inelastic scattering length 
exponent. In the limit of zero temperature, one expects p = It has been shown in 

experiments that the exponents are independent of which plateaux are involved, including 
the fractional ones. As we will see, our model has this feature. We will also argue that our 
model predicts the value: 



V = 20/9 



:i.3) 



There has been a certain prejudice for the value u = 7/3 = 21/9, though the only 



analytical computation supporting this value p9[ is based on the percolation picture of the 
network modelQ, and furthermore has been criticized as perhaps not corresponding to the 
right physics [0. Experimental measurements of z/ = 2.4±0.1 and u = 2.3±0.1 were reported 
in [30| and |31| respectively. If we are allowed to assume fi = 1/z/, then our result is rather 
close the most recent measurement of /i = .45 ±.05 reported in The numerical estimate 
with the smallest error, u = 2.35 ± .03 is due to Huckestein and includes data from both 
random Landau matrix approach to an Anderson tight-binding model in a magnetic field 
and from the network model A variant of the network model ||33| gives u = 2 A3 ± .18. 
Our value is closer to the result of Ando u = 2.2 ± .1 obtained from transmission through 

Indeed if the result 



a disordered system Indeed if the result ( p..3|) is actually correct then the consistently 
higher values for u obtained numerically for the network/ Anderson models suggests that our 
model is in a different universality class. Certainly it is clear from the very definitions of 
the models that they are not simply equivalent. 



'^v = 7/3 = Up + 1, where Up = 4/3 is the percolation exponent. 

^The data from the random Landau matrices and from the network model were statistically 
indistinguishable. (B. Huckestein, private communication.) 
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II. THE CRITICAL MODEL AND ITS CONFORMAL PROPERTIES 



A. Definition of the Model 



We consider fermions in two spacial dimensions with coordinates x, y, which we often 
denote simply as x. The constants e = c = h = 1. The second quantized action one needs 
to study is 

52+1 = j dtd^x {idt -H)^ (2.1) 

where H is the hamihonian. For the purpose of studying disorder in the x — y plane arising 
from a disordered potential in H, it is convenient to work with the time-Fourier transformed 
Green functions. Thus we expand 

^{x,t)= f ^ e-'^' ^> e{x) (2.2) 

and the action becomes 

^2+1 = Jde J dPx ^\{x){e - H)^,{x) (2.3) 

Since the functional integral is defined by e**^, for a fixed energy e one needs to study a 
euclidean functional integral defined by , where 

S = ij d^x "^^{H-e)^ (2.4) 

and we have dropped the subscript e on 

Our model consists of a two component fermion 

(2.5) 



with the hermitian Dirac hamiltonian 



H = {-id^ - A^) + ^ {-idy - Ay) ay + V{x, y) (2.6) 

This kind of Dirac hamiltonian has been considered before in connection with Quantum Hall 
transitions [|ll|] |jlO[, but here the meaning is rather different. Our theory is interpreted as 
describing the extended states in a bulk that is free of impurities. For spinless electrons it 
should be a scale invariant theory with c = 1 that is rotationally invariant, and the above 
hamiltonian is essentially the unique one with these properties. Wen's description of the 
edge states may be useful in making this connection more precise [Q. Defining the complex 
coordinates 

z = ^{x + iy), z = ^{x-iy), (2.7) 



and the gauge fields A^ = {A^ — iAy) / \/2, A^ = {A^ + iAy) / \f2, upon rescaling \1/ ^ v^27r 
one obtains 
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s 



(fx 



iPlid, - iA,)4j2 + 4^l{d^ - lA-)^^ - i{y + e) (^/'IV^s + V^iVi 



(2i 



The V and e terms give the fermions a mass and thus break the conformal invariance. For 
the remainder of this section we set = e = and focus on the conformal field theory 
coupled to the gauge field. 

The fermions can be bosonized with a single scalar field satisfying = 0. The 

scalar field separates into left and right moving parts 



(2.9) 



and the fermions have the bosonized expressions 

= e"*"^^, 4 = e*"^^, ^2 = e^"^^, = e'^'^^ (2.10) 

The current coupled to the gauge field then has the following bosonized expressions 



—iblihi = —dz 

277 



(2.11) 



chef) 



27r"'" 2ti 
Introducing the completely anti-symmetric tensor 



-xy 



—e 



yxi 



the current takes the form 



1 

2^ 



(2.12) 



(2.13) 



The current is thus a topological current which is identically conserved = by virtue 
of the anti-symmetry of 6^,^. The bosonized action takes the form 



S 



d X 



Stt 



27r 



(2.14) 



where d^df, = 8^3^ + d^d^. 

In order to clarify some of our subsequent arguments, we introduce a coupling g by 
rescaling the current in the j^A^ coupling to 



2'Kg 



-e„,ud, 



One justification for this is to consider adding a current- current interaction 

^ J d^xjf,]^ = k J (fx^{d^(t)f 



(2.15) 



(2.16) 



Such an interaction would arise for instance upon averaging over a gaussian disordered 
component of the gauge field A^, where k is proportional to the variance. This interaction 
merely re-scales the kinetic term. Redefining 0^0/ yjl + k leads to the current ( p.l5|) with 
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g"^ = 1 + k. If one prefers, g can be viewed as a kind of bookkeeping device that will be 
convenient in the sequel. Our original model corresponds to g = 1. 

Whereas the fermionic description has manifest gauge symmetry, the bosonic description 
as it stands does not. However, the gauge invariance can be restored by adding an term: 



S 



27ig 



e^udy(t)A^ 



2ng 



2 Afi^fi 



(2.17) 



To see this, introduce the dual field defined as 



In terms of left and right movers defined in Eq. (|2T 
constant) 



(2.18) 

the above equation implies (up to a 

(2.19) 



Using the relation [d^cf))" 
invariance 



-{dfj_<f)Y one can verify that the action has the local gauge 



0(a;) (j)(x) + -X{x), 
9 

Since ipi = exp(— i(0+0)/2), ip2 = exp(z( 
is 



A^{x) ^ A^{x) + d^X{x) (2.20) 
-0)/2), the gauge transformation on the fermions 



1,2 



-i\ 



1,2 



(2.21) 



for = 1, as it should be. The addition of the term is analogous to the Green-Schwarz 



mechanism for anomaly cancelation in low energy effective string theory [35 



For application to the Quantum Hall effect, the gauge field should incorporate some 
magnetic flux $ through the x — y plane. The flux is $ = / (PxV ^ A = § dx^A^. A 
non-zero flux can be obtained with a pure gauge, A^ = d^x^ loiig X is allowed to be 
discontinuous across a cut in the x — y plane extending from the origin to infinity. Namely, 
the flux is non-zero if x h^-s winding modes 



dad^X = = 27r) - x(0). 



T2 



(2.22) 



For Afj, = d^x the action takes the form 



S 



d^x 



2ng 



(2.23) 



The coupling between x ^i-nd is a topological term, i.e. a total derivative. A model of 
two independent scalar fields with the same topological coupling has been studied before as 
a c = 2 conformal field theory [|l^ [jl8|. Our model differs in a significant way: due to the 
structure of the couplings the gauge invariance Eq. ( p.20| ) allows x to be gauged away up 
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to the effects of its discontinuity across the cut. Thus we expect our model to share some 



of the features found in |17] |18], but in a c = 1 conformal field theory. 

As we now describe, a further constraint on the zero modes of x allows us to incorporate 
an essential feature of the 2 + 1 dimensional Quantum Hall dynamics. Electric current 
conservation implies the continuity equation dtp + d^J^ = where p is the charge density. 
When the conductivity axx = we may write = axy^^vE,,. Inserting this into the 
continuity equation and using the Maxwell equation V ^ E = —dtB, one obtains dt{p — 
a^yB) = 0, where B is the magnetic field perpendicular to the x — y plane. Integrating over 
space: 

I (g - axy^) = (2.24) 

where q = J d'^xp is the electric charge and $ = / d'^xB the flux. 

Let C denote a circular contour of arbitrary radius surrounding the origin. If the charge 
inside this circle is zero at t = — oo, then at time t. 



q{t) = - dti dx^e^^J, (2.25) 

J-oo JC 

For t oo, using the time dependence in Eq. (pl2|), 

dtJ^,{x,t) = j dej^{x) (2.26) 



where is the current in Eq. (|2.11 ). Since the flux is $ = §Qdx^d^x-, fo^' fixed e the 
equation ( |2.24| ) leads to 



dx^e^yjy = (Jxy £ dx^df^x (2.27) 



We take the point of view that the proportionality expressed in Eq. (|2.27| ) captures an 
essential feature of the Quantum Hall dynamics. Since this is real time dynamics in 2 + 1 
dimensions, it cannot follow from our action S, but must be put in by hand. We will express 
this proportionality in terms of a fundamental parameter 6: 



e f 



2ng 



j( dx^d^cj) = £ dx^d^x (2-28) 



where we have used Eq. ( |2.15| ). Though Eq. ( ^^.281 ) follows from Eq. ( ^^.27] ) with the 



identification 6 = l/a^y, we only expect this to be valid when = 0, and so it is incorrect 
to make this identification at this stage. Rather, in the next section we will compute the 
conductivities axx, <^xy in terms of the parameters g, 9. 

To summarize, our conformal model is defined by the action Eq. ( p.23| ) with the con- 
straint Eq. ( p.28| ) on the zero modes. 
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B. The Spectrum of Conformal Fields Depends on 9 



The action possesses the gauge invariance 



9 



(2.29) 



which follows from Eq. ( 2.20 ). Thus, x can be gauged away up to its effects on the zero 
modes. One consequence of this is that the spectrum of allowed fields is modified in a way 
that depends on 6. Consider the conformal fields exp(iai02. + ^«r0r)- The topological 
charge separates into left and right pieces: 



Q = - £ dx^e^J^ = Ql-Qr 



where 



Qi 



1 



dz dz 



Qi 



^ I dz ch(f)R 
Ittq J 



2'Kg J 2'Kg 
The exponential fields are characterized by their charges ul^r. 



(2.30) 



(2.31) 



(2.32) 



To determine the spectrum of allowed aL,R we conformally map the theory onto the 
cylinder by letting 



w = t + ia 



(2.33) 



The anti-symmetric tensor is e^t = —^ta = 1- The coordinate a is along the circumference 
of the cylinder and takes values < cr < 27r, whereas t runs along the length of the cylinder, 
— cxo < t < oo. In order to carry out canonical quantization, let us for the moment further 
rotate to Minkowski space t — > it. The field cj) can be expanded as 



1 



(j){t, a) = (t)o{t) + gQa + iJ2~ ('^"^ 



n 



in(t+cT) _|_ - ^-in{t-a) 



(2.34) 



The lagrangian for the zero modes is 



L 



2tt 



da 



in 



1 

2^' 



where $ is again the magnetic flux. 
The momentum conjugate to 0o is 



27rg 



idaX 



(2.35) 



1$ 



27rg 



(2.36) 
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Substituting (poit) = (f)o + (l>ot into Eq. (|2.34| ) and using Eq. ( p.36|) one obtains the zero 
mode contribution to (p. After mapping back to the z coordinate one finds 



4>{z,z) = 00 - iQi^ogz - iQRlogz+ ... (2.37) 

where 

Ql={po + ,^^]+'-^ (2.38) 
Qr = (po 

We first impose that the charge Q is an integer n. Then, since [^o^Po] = h mutual locahty 
of the exponential fields when ^ = requires po = m/g where m is some other integer. Finally 
the constraint Eq. ( p. 281 ) imposes $ = 6Q. Thus the spectrum of exponential fields consists 
of the operators 



* 


gQ 


27rg) 


2 




gQ 


27rgJ 


2 



On,m = exp (iaL0L + iaRCpR) (2.39) 



where 



if 9 \ gn 

"^ = ^r+2^"j-Y 

The integers n, m are electric and magnetic charges. 

The structure of the fields On,m can be understood as simply arising from the gauge 
transformation Eq. ( p.20| ). Namely, using Eq. (p.28| ) to identify x = ^0/2vr5f, the gauge 
transformation reads 

^^0+^0, 0^0 (2.41) 



Let us express 



aL0L + = a0 + a0 (2.42) 



where a = {ai + aRj/^, a = (a^ — aR)/2. Then, indeed one can verify that the 9 dependence 
in Eq. (g^ follows from the shift Eq. (g^: 

On,m{<P, 0; 0) = On,m (^<P, + 9 = 0^ (2.43) 

In summary, our model is a free massless scalar (p supplemented by the transformation 
Eq. ( 2.41 ). Any correlation function involving the fields (f),(f) will be computed by first 
performing the transformation ( |2.41| ) and then using the identifications ( p.9| )( prT9D and the 
two-point functions: 
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(0L(;2)0i(O)) = - log 2, (0fi(^)0R(O)) = -log^ (2.44) 
The identifications ( p.9|)(|2.19|) together with the above equation imply (when 9 = 0): 

(0(^,^)0(0)) = {^{z,z)m) = -^og{zz) (2.45) 
(0(z,^)0(O)) = {^iz,z)m) = -log ['^ 

Some 6 dependent correlation functions will be computed in the sequel. 

C. The Partition Function has an SL{2, Z) Symmetry 

The partition function on the torus possesses an SL{2, Z) modular symmetry acting on 
the coupling constants g, 9. Demonstrating this in our (c = 1) conformal field theory closely 
parallels the discussion in |TP[ ||TB[ for 2-boson (c = 2) theories with topological term. 

A torus is obtained by imposing periodic boundary conditions on the cylinder described 
above. Let the length of the cylinder in the t direction be /. Since the hamiltonian on the 
cylinder is H = Lq + Lq — c/12, where Lq, Lq are the zero modes of the Virasoro algebra, 
one has 

Z = Tr e-'(-^«+^0"^/i2) (2.46) 

The trace is over the Virasoro highest weight representations corresponding to the fields 
On,m- The zero mode contribution to Lq + Lq corresponds to the (anomalous) conformal 
scaling dimension of On,m- Using 

^g«ai</'i(z)g-jaL0L(O)^ _ ^giafl</'fl(z) g-«aH</>H(0) ^ _ (2.47) 

one finds 

dn,m{9^ 0) = 2 + = ^ h + ^ ) + V ^ 

The partition function is then 



Z{g, e) = ^Y. Hdn,m{9, 0)) (2.49) 

\v\ 



|2 



where the Dedekind //-function comes from the non-zero modes an-, an- (See 

Z is obviously invariant under 6 ^ 6 + 27r since this just shifts the integer m. It is also 
easy to show that 

dn,mi9,0) = d„ig',e') (2.50) 

with 

0'- ' 



(gyA + {e/27T) 
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Thus, Z is also invariant under [g, 6) — > {g' , 6'). 

Introduce a modular parameter r (not to be confused with the geometrical modular 
parameter of the torus Ttoms = il/2n) 

T = - + z'- (2.52) 
Zn I 

Then the two above symmetries correspond to 

T: + S : r ^ -1/r (2.53) 

These two transformations generate the group SL{2, Z), whose elements F transform r — *• 
r(r) = (ar + b)/{cT + d) where a, 6, c, ci are integers satisfying ad — he = 1. The partition 
function has the full symmetry Z{T{t)) = Z{t). 



III. CONFORMAL CONDUCTIVITY 

The conductivity tensor a^y has the following properties: axx = o'to, CFxy = —(^yx- In 
terms of the complex coordinate z, this implies 

(Tzz =<7zz = (3.1) 

The non-zero components are 

(^zz ^xx ~l~ '^'^xyi (^z'z (^xx ~l~ i(^yx is^'zz) (3-2) 

Thus it is natural to think of the conductivity as the single complex parameter a-^z- We will 
also use the covariant description 

(^xx 2*^/^/^' ^xy '^^^u^^fiu ('^•"^) 

The conductivity is usually expressed in terms of retarded and advanced one-particle 
Green functions, which can be computed from the action ( p.8| ) by including a small positive 
(negative) imaginary part to e for the retarded (advanced) Green function. Studying a 
localization/delocalization phase transition amounts to finding a renormalization group fixed 
point of the theory where the theory is conformally invariant. Our aim in this section is to 
understand the properties of the conductance at the critical point, thus we set the potential 
V to zero. The e-terms also break conformal invariance. In what follows we will derive some 
simple expressions representing conformal contributions to the conductivity when V = e = 0. 
In a sense these "conformal conductivities" represent the conductivity at a possible fixed 
point of a model. Later in the paper we will restore an impurity potential as a perturbation. 

We start from the Kubo formula in 2 + 1 dimensions. Throughout this section t is the 
real time in the 2 + 1 dimensional world. We work at finite temperature as a computational 
tool, taking the zero temperature limit at the end. The AC conductivity is given by 

a^,{u) = -U^,{u) (3.4) 

UJ 



12 



where n(a;) is an analytic continuation to real time of the euclidean (E) Matsubara quantity: 

U^,iuj) = U'-fJ{tu^u + zv) (3.5) 
where rj is small and positive, and 

U'-fJituj) = -J d^x jyTe"'^{J^{x,T)J,{0)) (3.6) 

Here, r = it is euclidean time, f3 is the inverse temperature, and the correlation is at finite 
temperature. 

In the Matsubara formulation we work with the action 

S= 1^ dr [ d^x ip\dr + H)^ (3.7) 



and expand the fields as follows: 

1|J{x,T)=Y,e-"'^^p^u{x) (3.8) 

V 

where z/ = 27r(r2 + l/2)//3 with n an integer. The action then takes the form 

S = Y.(3 fd^x i^l{x){H - iv)'4^,,{x) (3.9) 

For notational simplicity, we carry out the computation without displaying the spacial 
tensorial properties, restoring them at the end. The current is a fermion bilinear: 

J(x,r) = J2e'(-'^^)- ^L(x)^..(x) (3.10) 

Inserting this into Eq. (p.6|) one obtains 

U(^\^u;) = z(e-^ - 1) E / d'x^, (v;L(x)^..(x) J(0)) (3.11) 

There are well-known techniques for carrying out the sums over Matsubara frequencies 
(see e.g. [^). The sum over z/' may be performed by considering the contour integral 

where 

""^^^^ = -^PT^ ^3.13) 

and the contour of integration is a circle of radius R a.s R oo. In this way one picks up 
the poles in np{z) at z = i7i{2n + l)//3 with residue — 1//9. Assuming the only other pole is 
ai z = iu — iuj, one obtains 
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U'-^\zuj) = -(3e^-^ I rf^x^ (v;l_(x)^..(x)J(0)) (3.14) 
The summation over v is now studied by considering the contour integral 

/ 1^ ""^^'^ (C(x)V;,(x)J(0)) (3.15) 

The integrand is known to have branch cuts aX z = e + iuo and z = where e is real. Thus 
the contour must be chosen to run above and below the branch cuts along z = e + iuj ±i5 
and z = e ±i5. The other parts of the contour are along a circle at oo, thus the integral 
Eq.( |3.15D is a sum over three closed contours and all poles in np are picked up. The sum 



over u is then expressed as the sum of four real integrations over e coming from above and 
below the branch cuts. 

What is normally done is to factorize the correlation function in Eq. ( p.l4| ) into a 
product of two 1-particle Green functions; the result is an expression involving a product 
of the difference between retarded and advanced Green functions, i.e. the spectral density^. 
The retarded (advanced) Green functions correspond to taking 77 = 0+ ( 0~) in Eq. ( p. 51) . 
This would lead one to believe that when 77 = the conductivities are zero. However as 



explained in ||3^ a non-zero density of states requires a vacuum expectation value which 
survives when ?7 = 0, and an analogy has been made with spontaneous symmetry breaking^. 
Our aim is to extract the conformal (critical) contribution to the conductivity, i.e. the part 
that survives as e, 77 go to zero, so we follow a different procedure than the usual one. We are 
interested in the DC conductivity at = 0. As u; tends to zero, in particular when u = 26, 
the integrations along z = e + iu — i6 and z = e + i6 coalesce, and there is effectively only 
one branch cut at z = e. Integrations above and below this cut give: 

n(^)(2cu) = /^'e'"^ [d'x [ ^ {npie + tuj) - npie)) (J,(x) J,(0)) (3.16) 
J J 2m 



with Je = tpltpe where e has the same meaning as in Eq. (|2.3| ). We can now make the analytic 
continuation ( p.5|) and take the lj ^ limit. Using dsnp^e) = —S{e) at zero temperature, 
the integration over e sets e = 0. Recall e = corresponds to the conformal limit. Finally, 
rescaling J^^q[x) — > j{x)/l3, the zero temperature limit may be taken. Our final result is 
the simple formula 



-J (fix {j,{x)jM) (3.17) 



We denote the above quantity as a'^^, to emphasize that it is the critical contribution to a^^. 
(The disappearance of the l/27r in Eq. (|3.16|) comes from the extra 2tt in Eq. ( p.8|) .) The 
same rescaling of the currents hy 1/(3 removes the (3 in Eq. ( p.9|) . Thus in the formula Eq. 
( p.lTp , the currents are those in Eq. ( |2.11| ) and the correlation function is computed with 



^For a discussion in the context of disordered electrons, see [37|. 
^See also section IVC. 
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respect to the action ( p.23| ) with e = 0. The conformal currents j'^ have dimension 1, so 
that a'^^ is dimensionless as it should be. 

We now evaluate cr^j, in our model. The current is = id^(j)/27ig. We evaluate the 
current-current correlator in the presence of the coupling to x by simply performing the 
gauge transformation Eq. (|2.41|) : 

d'x Id, U{x) + —Mx)\ d, {m + —,m\ ) (3.18) 



We take the geometry to be a disk of radius r with r going to oo, and we don't impose 
any specific boundary condition at infinity. The integrals in Eq. (|3.18|) can be moved to the 
boundary of the disk using the divergence and Stokes theorems 

j cfxdf,f, ^ ~ / ^^f^^i^^f"^ j ^m;^ ^t^f" = / ^^m/m (3-19) 
Consider for instance the contribution: 

j2 



d'x (9^0(x)a^0(O)) = -j>^ dx, e^, (0(x)9,0(O)) (3.20) 
where the contour C is the boundary of the disk. Using Eq. ( p.45|) one finds^ 

/A(W.m) = 2.(/|ii + /i|i) =4. (3.21) 



Similar reasoning gives 

j d^x {d,4>{x)d,m) = 4vr, j d^x {d,4>{x)d,m) = (3.22) 

We also need 

' d^x e,, {d,4){x)d^4){Q)) = I dx,{4>{x)d,^{0)) = 4m (3.23) 



ic 

The same result holds with and interchanged. Putting this all together we obtain 

1 r /„, o\2 
ng^ 
i 26 



<. = ^(l + (^/V)") (3.24) 



a. 



27r7r^^ 



^The integral over d'z originally has opposite sign due to Ezz = —^zz, but for the contour C, the 
sense of integration over z is reversed in comparison to z; the integrals in Eq. ( p.21| ) are both usual 
Cauchy integrals. 
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The overall 1/2tt is expected from the fact that cr^y has units of /h = 1/2tt. 

That a'^y is imaginary for real 9 can be traced back to the "i" in the coupling to the gauge 
field in Eq. ( p.l7| ). We now give arguments supporting the idea that for the computation 



done in this section, we should analytically continue from euclidean to Minkowski space. The 
reason has to do with the difference between the hermiticity properties in first quantization 
verses second. In first quantization, the hamiltonian if is a hermitian operator H"^ = H. For 
our particular hamiltonian ( ^.61 ) this follows from {—18^^ = —id^ and Al y = Ax.y. However, 
when working with the action 5* as a conformal field theory, the natural reality properties 
of the functional integral are different from this. Note for instance that the bosonization 
( p.lOp is at odds with the hermiticity properties of il>i^2- Let us focus on the current 
in Eq. ( [2.11| ). In the first quantized viewpoint, (V^i,2)^ = "^1,2 implies {jz^ = jz, which 
requires (9^0)^ = 3^0. The latter is not the usual hermiticity assumed in conformal field 
theory which is tied to the functional integral. In conformal field theory one considers rather 
{dz4'V — ^z<P- This suggests that to perform the computations in this section meaningfully, 
one should analytically continue to Minkowski space. The complex analytic structure of the 
conductivity tensor as displayed in Eq. ( p.2|) also allows the interpretation of the analytic 
continuation of continuation from euclidean to Minkowski space for the coordinates 

X, y. From the point of view of first quantization, this seems to correspond having an 
imaginary vector potential. We point out that the importance of imaginary vector potentials 
for delocalization transitions has recently been recognized in both electronic and biological 



systems [3£] [40 [41|. For more discussion on this rather delicate point, see [14 



We define now the critical condition 

= ±9/2% (3.25) 

There are two interpretations of this condition depending on the prescription taken to make 
axy real, as we now describe. 

We can perform the analytic continuation to Minkowski space by letting 9 —i9. The 
critical condition then leads to a' = o"^ with 

= 0, = 1/9 (3.26) 

Note that this is precisely the identification we made in arriving at the zero mode constraint 
Eq. (p.28|) . In terms of r, the critical condition leads to 

r^ = |-(l±z) (3.27) 

For our original fermion model with g = 1, the critical condition gives 9 = ivr, so that 
r'^ = {l + i)/2, {—l+i)/2. These are known to be the non-trivial fixed points of the SL(2, Z), 
one being related to the other by T. Our interpretation of the conformal field theory with 
the critical condition is that it represents a pure, ideal system consisting only of the extended 
states, which are the extension of the edge states into the bulk. This pure system is known 
to have the same conductance properties as a physical plateau in the presence of impurities 
0, so Eq. ( 3.26| ) leads us to this interpretation. We emphasize that we have not added any 



impurities yet, so that even though Eq. ( |3.26| ) has the same properties as a plateau, there 
are no real plateaux yet in our model since there are no localized states. We cannot speak 
of the critical values of a^x, <^xy at a transition until impurities are added. 
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The other possibihty for making axy real is to analytically continue y ^ iy, keeping 9 
real as before, under which a^y — > 9/Ti'^g'^, and a^x remains unchanged. Here the critical 



condition still leads to r*^, and 



al^ = ±al, = 1/9 (3.28) 



For this prescription of analytic continuation, since ^ 0, the interpretation would be 
that we are building into the model some features of the transition with impurities, i.e. the 
property that ^ 0. Since we will next add impurities as a perturbation, the previous 
interpretation is preferable. 

The SL{2, W.) symmetry ( |2.52| ) does not have a simple action on the conductivities since 
the latter are not a modular transformation of r. We can however point out the following. 



We are mainly interested in couplings satisfying the critical condition ( |3.25| ). Let us choose 
the positive sign in ( p.27|) , and construct a modular parameter out of a^x, CTxy' 



^ = 2ma,^ (3.29) 

Then one has: 

? = 5T"(r'=), for 9 = -27m (3.30) 

So, for this restriction of the parameters, modular transformations of r induce modular 
transformations of and this should lead to certain features of the phase diagram. 



IV. ADDING A CIRCULAR DEFECT OF IMPURITIES 

A. The Impurity Potential 

We now introduce a potential V representing some impurities in the system. For the 
reasons explained in the Introduction we do not chose a random potential. We will take the 
potential to be of the form 

V{x,y) = VoS{r-ro) (4.1) 

where = x"^ + y^ and Tq is some arbitrary radius. This corresponds to a circular defect 
line of impurities at r = tq. As before, the sample is a disk of radius r — oo. This poten- 
tial certainly incorporates a fixed realization of impurities. As we will see, the important 
advantage of our choice of potential is that the critical exponents do not depend on Vq. 
The potential ( [4.1D leads to a term in the action Eq. (|2.8| ) 

72 

Sy = -tVo l^^ir- ro) Oy(x, y) (4.2) 



After bosonization, the model we will study has the action Eq. ( |2.23| ) with the additional 
perturbation Sy, where Oy = cos (p. 
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B. The Defect Theory can be Mapped into a Boundary Field Theory 



Let us map the theory to the euchdean cyhnder as in section II, where the coordinates on 
the cylinder t, a are defined in Eq. (|2.33|) . Since = +?/^ = 22:2 = 2e^\ the defect IS now 
a circle along the circumference of the cylinder (the a-direction) at to satisfying tq = \/2e^° . 
Without loss of generality we let to = 0. 

The theory can be solved by folding it onto a boundary field theory. This folding proce- 
dure has been previously applied to other physical problems |^ |^ We first set 
Vo = and fold the conformal field theory. In the defect version of the problem any field 
can be separated into its pieces on either side of the defect: 

0(t, a) = e{t)<P^+\t, a) + e{-t)(j)^~\t, a) (4.3) 

where 9{t) is the step function, 9{t) = 1 for t > 0, zero otherwise. Let (t)Lji denote the L — R 
components of the scalar field on either side of the defect: 

0W = 0W(^ + ,^) + ^(±)(^_,^) (4.4) 

From these we define the following fields for t > only: 

ip = ipL + ipR, LpL{t,a) = (j)-^\t,a), V9ij(t,a) = 0if^(-t,o-) (4.5) 

^' = ^'l + ^'w ^'lit^ ^) = 4>R\-t, a), <^'^(t, a) = 4>'^^\t, a) 

As defined, the fields Lpi.Lp'j^ {Lp^.ip'j^ are functions oi t + ia (t — zcr), hence their L/R 
designations. Since they are defined only for t > 0, they are fields in a theory with a 
boundary at t = 0. Define now the even/odd combinations: 

<^(^) = ,p + ^\ <^(°) = ^-v' (4.6) 

The action for the boundary field theory corresponding to the free hamiltonian can be 
written as 

^free = - dtj da- {{d.^^^'^f + {d,^^"^f) (4.7) 

The theory is not fully defined until the boundary conditions in the conformal field theory 
are specified. General properties of boundary conformal field theory were studied by Cardy 
47| . We first set the coupling to the gauge field to zero. In the defect description, when 



V = Q the appropriate boundary condition is dt(j) = since this corresponds to being 
continuous across the defect. Since (pi, 4>R are analytic functions of w = t + ia, w = t — ia 
this implies that on the defect (pL + 0r is a constant, which we take to be zero. Letting 
= (0*-^-* + (j)^~^)/2 on the defect, in terms of the boundary fields 0l + 0_r = reads 

+ <PV + 0« ^ + <PR^ = ^^'^ = 0, at t = (4.8) 

Consider next the impurity operator Oy- In the defect formulation it is appropriate to 
take COS0 at the defect to be the average of its values on either side of the defect: 
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COS0 = ^ (cOS0(+) + COS0(-)) = COS + 0^"V2] COS [(0^+) - 0^"V2] , = 0) 

(4.9) 

Translating this to the boundary description, on the boundary one has 

0(+) - = + - _ = - ^« - + ^'^ = (^(°) (4.10) 

where as before the dual field (^(°) is defined as d^i^^"^ = —ie^,^d^ip^°\ Thus on the boundary 
one has 

Ov = cos((^(°V2), {t = 0) (4.11) 

Since the field (f^^^ decouples from the boundary, we henceforth drop it. In order to 
match the normalization of the previous sections we scale out the factor of 1/2 in Eq. ( |4.7|) 
by now defining 

= (^(°Vv^ (4.12) 

The gauge field can now be restored; for = the x terms can also be folded as for 
the field. We finally obtain the action for the boundary theory: 

(4.13) 

where in the above equation x = The coupling g is related to g as follows 

g = V2g (4.14) 

Finally we need to impose a zero mode constraint as in Eq. ( |2.28|) . The original charge 
Q in the defect theory can be expressed as Q = / dadacj)^^'' /2Txg. At t = this is Q = 
/ dada{'^^^^ + (p^°^) /4:7rg. Dropping (^^'^^ and using Eq. ( ^.12| ) one has Q = J dad„(p/27rg. 
Thus we impose 

^ <j) dx^d^ip = j dx^d^x (4.15) 

where 6 is the same as in the unfolded theory. (Here x^ denotes t, a.) The couplings g, 6 
are real, which leads to real anomalous dimensions of operators, and the critical condition 
is (Ipl). 

The original fermion model is characterized by the values g = -\/2, 6 = n. When ^ = 0, 
g = \f2 is the self-dual point of the SL{2, Z) symmetry described in section II, and Vq 
corresponds to a (boundary) marginal perturbation of scaling dimension 1. 

Our model still possesses the gauge invariance ( |2.29| ) with g replaced by g. A 2-boson 



model with topological coupling and a similar boundary interaction was considered in 
again the important difference is the gauge invariance of our model. 



19 



The new feature of the boundary version in comparison to what we had in section II 
is the boundary condition at t = 0, which now depends on 9. Letting ip ^ ip + dip and 
requiring (55* = on the boundary leads to the equation of motion (when = 0): 



Imposing the zero mode constraint Eq. (|4.15| ) one finds 



(4.16) 



f/2 + 6/271) d^^L 



f/2 - 6/21,] d^^R {t = 0) 



(4.17) 

This is an interesting relation since it links the modular properties of the couplings 6 
with the spacial coordinates. Namely, going to Minkowski space and defining the light-cone 
coordinates v = t + a^v = t — a, the boundary condition reads 



{rd^ - rd^) if = 



(4.18) 



where r is the modular parameter ( ^.52| ) with g ^ g and r = r*. 

We will need the anomalous scaling dimension of the boundary operator Oy- We first 
gauge away x'l before the zero mode constraint leads to the transformation ( p.41[ ) with 
g ^ g so that 

(4.19) 



(4.20) 



Ov ^ cos ([(l + 6/7if) <^L - (l - 6/ng^) ifR 
Viewing a as the "time" the hamiltonian can be written as 

H={^ r dt idn,^L? + {d^<l>Rf) + ^ Ov{t = 0) 

VOTT JO / ZTT 



The theory can now be rewritten using only a left-moving field. Using arguments found in 
47| , the boundary condition Eq. ( 4.17 ) allows us to view ip^ as an analytic continuation of 



ipL- Using Eq. (|4.17| ) we make the substitution 

if /2 + 6/271) 
{gy2 - 6/2n) 

into Eq. ( [4.201 ). Defining a rescaled field 



(4.21) 



(1 - 6/7rg'^) 



one finds 



where 



H 



1 



d(y{d^(pL) 



2t\ 



■ cos 



V2a(pi 



(4.22) 



(4.23) 



(4.24) 
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Since the rescaled field has the two-point function {(Pl{z)i^l{0)) = — logs;, and 

^^iV2a^L{z)Q-iV2a^L{0)'^ _ J^/^jSa^ ^ (4-25) 

the anomalous scaling dimension is 

[Ov] = dim(Cy) = (4.26) 

C. Critical Exponents 

The model of the last section has a single energy scale M with units of mass set by 
Vq. Since the action is dimensionless, the dimension of Vq is 1 — [Oy]- The energy scale M 
defines a correlation length M = 1/C,c, thus 

oc (4.27) 

Let p{E) denote the density of states, so that p{E)dE represents a number of states per 
unit volume. It can be expressed as 

p(E) = ^Tr 5(H - E) = ^ lim Im Tr ( — ^ ) (4.28) 

where V is the 2d volume. This in turn can be expressed as a retarded Green function 

p{E) = - lim Im Ul{x)^s{x)) e = E + ir] (4.29) 

TT »7— >0+ \ ' 

In accordance with the discussion in section III, we define a critical density of states p'^ from 
the above formula with e = 0. As we argued in section III for a'^^, this quantity should 
represent the density of states near the critical point. Since both fields in Eq. ( [4.29| ) are at 
the same point x, p'^ is a one-point function of the operator ip'^il) = iplipi + ip24'2- From Eq. 
( ^72|) this operator is the operator Oy- 



oc {Ov) (4.30) 

From the scaling dimension of {Ov) we know that 

{Ov) oc {a-^""^^ (4.31) 

Next, recall that the number of states per unit volume in a Landau level is 3/271. Thus, 
p'^ should scale with 5, so that p'^ oc \B — B^]. Using Eqs. (|4.30|)(|4.31|) one obtains the 
relation ( p,.2p with 

u = l/[Ov] (4.32) 

Alternatively we can argue that since since p{E)dE has units of inverse volume, in two 
dimensions p has dimensions of energy. Letting p'^ oc E, one then finds 

e^oc (4.33) 
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where here Ec = 0. The connection between the two equations ( [4.33|) and ( [1.2|) is due to 



the fact that the Fermi energy is proportional B in a system with Landau levels. 

In order to satisfy the critical condition a'^^ = 0, from Eq. ( p.25| ) we impose O/ng'^ = ±1. 
Using Eq. ( [4.19|) one has O/irg"^ = ±1/2. Inserting this into Eq. ( [4.26|) we obtain the 
value u = 20/9 quoted in the Introduction. Note that this exponent only depends on the 



unconstrained, and this indicates that the exponent u is universal. 



ratio 0/g . This value of v holds for any g,6 with a'^^ = 0. In particular a'^y = 1/9 is left 



V. S-MATRICES DESIGNED FROM THE MODULAR PARAMETER r 
EXHIBIT A SERIES OF PLATEAUX 

In this section we propose an intriguing connection with the so-called staircase model, 
whose interesting properties were first witnessed by Al. Zamolodchikov [4£], and were con- 
sidered mysterious at the time. 

The conformal model of section II is essentially that of a single massless scalar field with 
coupling constants defining a modular parameter r and possessing an SL{2, Z) symmetry. 
The original staircase model was characterized by a bulk S-matrix for a single massive 
particle. Though it is the boundary version that is more appropriate to our problem, let us 
begin by showing how the structure of the S-matrix follows quite naturally from what we 
have done. The energy and momentum of a relativistic particle can be parameterized by a 
rapidity /5: 

£' = mcosh/3, P = msinh/3 (5.1) 

where m is the mass of the particle. The bulk S-matrix describing the 2-particle to 2-particle 
scattering must satisfy crossing symmetry and unitarity ||50[] : 

Si/3) = Sitn - (3), Si(3)Si-(3) = l (5.2) 

S-matrices satisfying the above functional equations are usually built out of products of the 
minimal factors: 

QfQ\ sinh/3-2sin7r7 

S{f3) = . . (5.3) 

smh p + t sm ir^y 

where 7 is a parameter related to the coupling constants of the theory. Let us attempt to 
relate the S-matrix to the physics of our problem by using the SL{2, Z) symmetry to relate 
7 to r. The S-matrix has the following symmetries: 

7-^7 + 2, 7^1-7 (5.4) 

It does not seem possible for the S-matrix to possess the full SL{2, Z) symmetry. The first 
symmetry in the above formula suggests the identification 7 = r, since it would correspond 
to the transformation T^, where T is defined in Eq. ( |2.53| ). Next consider the transformation 
7 — 1 — 7. The modular parameters r satisfying the critical condition Eq. ( |3.25 ) come in 



complex conjugate pairs Eq. ( |3.27| ). It is natural then that the S-matrix not distinguish 
between these two critical r's. Requiring then that 7 ^ 1 — 7 is equivalent to rt 
requires = 1, or equivalently ^ = vr, in which case 7 = = (1 ± Finally we perturb 



22 



away from the critical condition by deforming g away from 1 in r, but keeping 6 = tt. 
Because of the form of the critical condition, we can just as well view this deformation as 
a deformation of 6 away from vr. Thus we write this as — s> 60/2-11, and identify 7 as 
follows: 

1 9 

^ = r± = = TT, = ±60/211) = -±i^ (5.5) 



2 27r 



The S-matrix is then 



sinh P + i cosh 60/ 2 

It was discovered by Al. Zamolodchikov that the above S-matrix leads to a free energy 
with some remarkable properties. The model can be studied on a cylinder of length / and 
radius R. Imposing periodic boundary conditions in the /-direction, the free energy S{R) 
was computed starting from the S-matrix by means of the thermodynamic Bethe ansatz 
|5T| . It was found in that as one varies R the free energy S{R) goes through a series of 



plateaux, i.e. the free-energy as a function of R resembles a staircase. 

Let us try and be more specific now about the relation with our model. The boundary 
theory we obtained is defined by the action ( [4.13|) . Let us first decouple the gauge field, 



X = 0. Noting that {d^ipY = — (S^jC/?)^, letting ip ^ i(f) we obtain from ( [4.13| ): 

Vo 



00 



S = dt da 







^(9^0)2 + Acosh(^</.) 
an 



^ I da cosh(/30/2) (5.7) 
2111 J 



with (3 = \/2. Our model has zero bulk coupling A = 0. We have included a A term since as 
written the above action defines the boundary sinh-Gordon model. It is an integrable model 



that can be studied using the framework developed in The bulk S-matrix is known to 
be of the form Eq. ( |5.3| ) with 

7 = (5.8) 

First note that our model leads to /? = \/2 and 7 = 1/2. This is precisely the same as in 
Eq. ( ^.5|) when 6*0 = 0. It is therefore clear that our model ( [4. 131) is very closely related to 
the boundary sinh-Gordon model with = ^J2 with A = and 7 analytically continued as 
in (|5.5|) . This analytic continuation of 7 incorporates the topological coupling to the gauge 
field (9X) and amounts to continuing the coupling g^ /2 to the modular parameter r. In terms 
of /?, the above analytic continuation of 7 corresponds to a simple phase: 

A boundary version of the staircase model was studied in |^.[] Indeed it was found that 
the boundary entropy reveals a series of plateaux as a function of renormalization group 
scale. 



^This paper also includes a useful reproduction of the bulk results in 
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Keeping within the scope of this paper we refrain from going any more deeply into the 
integrabihty structures that allow a detailed study of this proposal, but will return to this 
in a future publication. 
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VI. DISCUSSION 



In summary, we have constructed a reasonably simple model which we argued describes 
the critical properties of Quantum Hall plateau to plateau transitions. It is a c = 1 con- 
formal field theory mainly characterized by a gauge symmetry and the SL{2, Z) symmetry 
which follows from it. Adding impurities singles out a particular operator Oy, whose scal- 
ing dimension, as computed in the conformal field theory, follows largely from the gauge 
invariance, and leads to the exponent u = 20/9. 

Experimental errors at this time are perhaps too large to distinguish between u = 2.3 and 
u = 2.2. However if the measurements continue to indicate values closer to z/ = 20/9, as in 



34 , this suggests that the transition is in a different universality class than the percolative 



class of the network/ Anderson models, assuming the numerical work on the latter is correct. 



Our result appears to be more consistent with the model used in Ando's simulation p4 
which gave u = 2.2 ± .1. 

Since the gauge symmetry is an important feature of our model it can perhaps be viewed 
as a simplified version of the Yang-Mills theories in higher dimensions which exhibit the 
S'L(2, Z) electric/magnetic duality pl| . 

The boundary staircase model described in section V is a promising candidate for a 
model that exhibits a series of plateaux transitions. 
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VIII. ADDENDUM 



After the first version of this paper appeared, we reported a study of the relevance of 
disorder in 
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There we gave evidence for two possible universality classes in the presence 
of disorder, in one of which all disorder is driven irrelevent by the presence of disorder in 
the gauge field. 
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